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1. Introduction 

Effective range function plays central role in the analysis and interpretation of two- 
particle low-energy scattering data. It is well known from the simplest case of s-wave 
scattering by a short range potential that it is more advantageous to study the effective 
range function k cot S than the phase shift S. The former is an even function of the 
momentum k and close to threshold can be expanded in a power series of k 2 , known 
as the effective range expansion. The first two coefficients in this expansion: the 
scattering length and the effective range, respectively, carrying important information 
about the underlying interaction, are regarded as fundamental parameters in low- 
energy scattering phenomenology. The inclusion of the long range Coulomb interaction 
is a non-trivial complication because the latter changes dramatically the singularity 
structure of the scattering matrix. The Coulomb corrections to the effective range 
parameters are not only model dependent but can be also quite sizable. Especially 
sensitive is the scattering length where, strictly speaking, there is even no upper limit 
to the size of the Coulomb correction, i.e. the scattering length that is finite in absence 
of Coulomb interaction, could become infinite when Coulomb interaction is turned on. 
Perhaps the best illustration that in real life Coulomb correction can be quite large 
is provided by the nucleon-nucleon scattering. With exact isospin invariance the nn 
and pp 1 Sq scattering lengths in absence of Coulomb interaction should be equal but 
experimentally |, |]: a nn = -18.6 ± 0.5 fm and a pp = -7.828 ± 0.008 fm. It has 
been confirmed by model calculations that most of the difference in these two values 
may be attributed to the Coulomb interaction. 

The analytic continuation of the effective range function to negative energies 
presents considerable interest as it provides simple means to locate the near-threshold 
singularities of the scattering matrix, that might be associated with shallow bound 
states. For conducting a direct calculation of the effective range function below 
threshold, the standard Coulomb wave functions for negative energies are required, 
and, efficient algorithms for calculating them are currently available in the literature 
j|, |j. Nevertheless, this approach leads to difficulties as the standard Coulomb wave 
functions are ill chosen for the purpose and for attractive Coulomb potential they are 
singular at the Coulomb bound state energies. Therefore, the usual procedure employs 
the effective range expansion for effecting the extrapolation to negative energies. 

In this work we wish to present a direct method of calculating the effective range 
function below threshold based on the solution of the Lippmann-Schwinger equation. 
In this equation the Coulomb potential does not appear explicitly but has been 
accounted for exactly by introducing Coulomb-modified Green's function and Coulomb 
distorted ingoing wave. By selecting a specific Coulomb Green's function, intimately 
connected with the K-matrix, we have been able to devise a calculational scheme that 
is free from the Coulomb singularities. Another advantage of the proposed method is 
that it provides for a simple solution of the Lippmann-Schwinger equation. This is a 
consequence of the fact that for negative energies the Coulomb Green's function has a 
separable representation of the form of a Sturmian expansion so that the kernel of the 
integral equation becomes degenerate. The task of solving the Lippmann-Schwinger 
equation is thereby reduced to that of solving a system of linear algebraic equations. 
The method can be applied either for calculating the scattering matrix, or the K- 
matrix below threshold. In particular, this approach can be used for locating the 
poles of the scattering matrix what is equivalent to solving the bound state problem. 
The K-matrix yields the effective range function and allows to extract the low-energy 
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scattering parameters. 

The plan of the presentation is as follows. In the next section we outline 
the Lippmann-Schwinger equation formalism necessary for calculating the scattering 
matrix as well as the K-matrix. In section 3, we employ the Sturmian expansion of 
the Coulomb Green's function for converting the underlying integral equations into a 
system of linear algebraic equations. Finally, in the last section we discuss possible 
applications and examine numerically the performance of the presented method. 



2. Effective range function 

Consider a two-body, two-potential quantum mechanical scattering problem with 
spherical symmetry where the interaction is a superposition of the Coulomb potential 
plus a short ranged nuclear potential V(r). It will be assumed in the following that 
the nuclear potential for large separations r has at least exponential fall-off ~ e~ Ar 
with A > being the inverse range parameter of the nuclear force. For small r, 
we assume that the behaviour of V(r) is not worse than r~ 2 . The Coulomb potential 
corresponding to point-like charges is taken as a Z/r where Z is dimensionless strength 
parameter Z < 0(Z > 0) for attraction (repulsion) and a is the fine structure constant. 
The effective range function is defined as the inverse of the K-matrix Ke(k) and is 
related to the scattering matrix by means of a generalized Heitler formula 

K e {k) n{k) 

where I denotes the orbital momentum and the independent complex variable k on 
the positive real axis becomes the physical cm. momentum and the complex function 
ge(k) can be interpreted as a generalized barrier penetration factor. This function has 
been studied in several papers §, |, |, |, 0, ||, and, g t {k) = i k 2e+1 /[(2£ + l)!!] 2 for 
Z = 0, whereas for Z ^ its explicit form is 

gt(k) = {2 m Zf i+1 /[{2i + l)!] 2 fl (1 + m 2 Ay 2 ) 

x {± + in) + V(l - h?)] - (tt/2) cot(i^) - log(ki|Z|/£0} , (2) 

where rj = fiaZ/k, /i is the reduced mass of the system and iJj denotes the logarithmic 
derivative of the gamma function. We choose the units system where h = c = 1. In 
(Q) the complex function T^(fc) is the scattering matrix that will be obtained from the 
solution of the Lippmann-Schwinger equation incorporating the appropriate boundary 
condition 

/>oo 

u t {k,r) = (j){k,r) + {r\G+{k)\r')V{r')ui{k,r')dr', (3) 
Jo 

where Ui(k,r) is the sought for wave function. In order to account for the Coulomb 
interaction, in (Q) the ingoing plane wave has been replaced by the regular Coulomb 
wave function ^(fc,r) and, respectively, the free resolvent by the Coulomb Green's 
function (r |G^"(fc)| r') satisfying the outgoing wave boundary condition. It is easily 
checked that acting on both sides of (||) with the Coulomb Hamiltonian, we retrieve the 
Schrodinger equation appropriate for the two-potential problem under consideration. 
Given the solution of (Q), the scattering matrix Ti(k) is obtained from the formula 

/•OO 

T t (k) = -2n Mk,r)V(r)u e (k,r)dr, (4) 
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where the regular Coulomb solution <fit(k, r) satisfying the following boundary 
condition at r = 

lim r~ e ~ 1 (pe(k,r) = 1, 

is well known [Io| and can be expressed analytically as 

&(fc,r) =r e+1 e ±ikr 1 F 1 (£+l±i V ,2e + 2;T2ikr), (5) 

where 1 F 1 (a,6;z) denotes the confluent hypergeometric function ]l2|| . Since the 
boundary condition is independent of fc, it is evident that 4>i(k,r) must be entire 
analytic function of fc 2 for every non-negative value of r and is real for real fc 2 . 

For physical momenta fc the function Ti(k) is related to the Coulomb modified 
nuclear phase shift Se 

1/T t (k) = Cui-qf k 2t {k cotSi -ifc) (6) 



where 



and 



Co{v) 2 = 2^/ [exp(2^) - 1] . (8) 

Although for physical momenta k the functions 1/r^k) and ge(k) are both complex 
but their imaginary parts are equal in magnitude and have opposite signs. Indeed, in 
this case from (||), one has 

gt{k) = Ctijif k 21 {2 M aZ[Re^(l + \rj) - log(|»j|)]/Co(»?) 2 + ifc} (9) 

so that using (^|) and (^|) in (H), the imaginary parts are cancelled and the resulting 
effective range function (Q) is real, as it should be. However, the effective range 
function is known to be real analytic function of fc 2 , so that it must be real also on 
the imaginary axis. For fc located on the positive imaginary axis the Green's function 
and <fii(k, r) are both real what implies that r^(fc) is necessarily real as well. Also the 
function ge(k) is real on the imaginary axis 

£ 

g t (k) = (2 m Zf^l\(2l + l)!] 2 [] (1 ™ 2 lf ) 

x {i[V(l + ^)+V(l-^]-(7r/2)cot(^7r)-log(|^|)}, (10) 

with f] = jiaZ/\k\ and this ensures that the ensuing effective range function is real. It 
should be noted that in ( fio| ) the poles in -0(1 — fj) ar e cancelled with those from the 
cotangent term and only the poles coming from ^(1 +fj) will survive, they correspond 
to the Coulomb bound states. 

To complete the above scheme for calculating the effective range function we need 
the explicit formula for the Coulomb Green's function. The latter involves a product of 
the regular Fi(rj,p) and the irregular Ge{r],p) Coulomb wave functions (with p — kr) 
defined in jfj), and takes the form 

{r\G+{k)\r') = -{2^/k)[G t {r 1 , P> )+iF l { n ,p > )] F t (i,, P< ), (11) 

where F e (r),p) = C e [rj) k l+1 <t) e (k,r). 

Actually, the K-matrix can be calculated directly, i.e. without reference to the 
scattering matrix T({k). To this end, a different Green's function that is real analytic 
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function of k 2 will be introduced. Such Green's function, denoted hereafter as Gf to 
emphasize its connection with the K-matrix, can be easily constructed and reads 

(r \Gf (k)\ r') = (r \G+(k)\ r') + 2ng t {k) &(fc,r) <j>z{k,r'). (12) 

Indeed, the above Green's function has all what it takes. For physical momenta the 
imaginary parts in the two terms on the right hand side of ( |l2| ) cancel each other and 
(r |G^"(fc)| r') is real. For k values located on the positive imaginary axis both terms 
are real and so is (r \ Gf {k) \ r'). Using spectral representation of (r \G^{k) \ r'), it 
is easy to verify that the Coulomb bound state poles, present in (r |G/(fc)| r') for 
attractive Coulomb potential, will be exactly cancelled with the corresponding poles 
that occur in gt{k) (cf. (|Io|)). As a result, the Green's function (r \Gf{k) \ r') given 
by (|l2|) is free from the Coulomb singularities and is a real analytic function of k 2 in 
any finite part of the complex k-plane for non-negative r and r' . The K-matrix will be 
obtained from the solution of a Lippmann-Schwinger equation that parallels equation 



w t {k,r) = <f>(k,r) + (r\Gf(k)\r')V(r')w e (k,r')dr', (13) 



o 

where Wi(k,r) is a new wave function that differs from ue(k,r) by a constant factor. 
Given the solution of (^3|), the K-matrix is calculated from the formula 

poo 

K e (k) = -2n cj> e (k,r)V(r)w e {k,r)dr, (14) 
Jo 

and the above function is meromorphic in any finite part of the k-plane. It is easy to 
check that setting we(k 1 r) = [1 — ge(k) Ki{k)]ui{k,r) in ( |l3| ) and using (||) - (||), we 
immediately retrieve the Heitler formula ([j]). This shows that the Green's function 
( p^| ) is unique for a given choice of ge(k). It should be recalled here that gi{k) is 
not free from ambiguities and can be determined only up to a polynomial in k 2 (cf. 
[jz], |^|) but the latter is normally set equal to zero. This has been also our choice in 
the present work. 

The Green's function (|l2|) can be also cast into an alternative form which is on a 
par with ( JTT| ) 

(r\Gf{k)\r') =-2vMk,r<)6e(k,r > ), (15) 

where the imposed asymptotic boundary condition has been made explicit. The 
irregular Coulomb wave function 9i(k, r) present in (|l5|) has been introduced in || M . 
The latter is a real analytic function of k 2 for non-negative r and can be expressed as 
a superposition of the standard Coulomb functions 

9 t {k, r) = Cfin) k l G t { n , p) + Reg e (k) cp t {k, r). 

Admittedly, the asymptotic behaviour (r — > oo) of this function is rather complicated 
but it is worth to pay such a price in view of the fact that the Green's function given 
by (|l5| ) is manifestly a real analytic function of k 2 , and, as such, can be systematically 
expanded in a Taylor series in powers of k 2 . The form ([l5|) is particularly useful for 
calculating the K-matrix at k = 0. Following (To), one has 

" [{21 + l)!/^ +1 ] v /^7/ 2 , +1 (2 v /^7) for Z > 

<M0,r) = { r e+1 for Z = (16) 

[(2£+l)!/^+ 1 ] v /^7 J 2l+l {2^pr) forZ<0 
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where (3 = 2[ia\Z\, and (J21+1, 121+1) ar e Bessel and modified Bessel functions, 
respectively [Q. The function 9i(0,r) is given in terms of the Neuman Y^+i and 
modified Neuman function K 2 e+i, respectively [[l2] 

( [P e /(2i + l)ly/j3? K2t +1 (2y/p?) for Z > 

0<(O, r) = i r~V(2^ + 1) for Z = (17) 

[ -tt/2 [/?V(2^+ l)!]v^37y 2 , +1 (2V^) for Z < 

Since the K-matrix is a real analytic function of k 2 the effective range function 
can be expanded in a Taylor series in powers of k 2 around threshold 

1 =±— + ±nk 2 + ---, (18) 



Ke(k) at 

where the expansion coefficients for £ = have their standard meaning, i.e. ao is the 
scattering length and ro is the effective range and the higher order terms coefficients, 
that have been suppressed in (|l8|), are known as shape parameters. Finally, the 
sign ambiguity in (JL8|), reflecting the coexistence in the literature of two different 
definitions of the scattering length, has been resolved by adhering to the generally 
accepted convention that (— ) sign applies for baryon-baryon scattering and (+) sign 
in all remaining cases. 

When the scattering matrix has a bound state pole at k = kg, then from ([!]) we 
obtain an equation for the binding energy B — kj 3 /2(i involving the effective range 
function continued below threshold 

For shallow bound states the effective range function can be continued below threshold 
by employing the expansion (|l8|). Formula ([[9]) correlates then experimentally 
measurable quantities: the low-energy scattering parameters and the binding energy 
in a model independent way. 



3. Sturmian expansion method 



Following the scheme devised in the preceding section, the scattering matrix and 
the K-matrix will be obtained by solving a dedicated Lippmann-Schwinger equation. 
Actually, the solution of the latter below threshold presents a quite simple task thanks 
to the Sturmian representation of the Coulomb Green's function. 

The Sturmians, or Sturm-Liouville functions |l3|, |l4|, [l^] Si u (b,r), where 6 is a 
free parameter (Re b > 0) are enumerated by the nodal quantum number v — 0, 1, 2, ... 
and result as solution of the differential equation 

' d 2 _ 1(1+1) 26(^ + 1 + 1) 
dr 2 r 2 r 



-b 1 



Siu(b,r) = 0, 



(20) 



with the following boundary conditions imposed, respectively, at zero, and at infinity 
]im S iv (b,r) (21) 

lim StJb t r) ~ e - br . (22) 

r — >oo 

The eigenfunctions are given in a closed form 



St„(b, r) = \v\j{2i + 1 + v)\} 1/2 (2br) 



- br Ll t+1 (2br), 



(23) 
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where L^ e+1 denotes Laguerre polynomial. The functions (|2^) are orthogonal and 
form a complete set with respect to the weight function 1/r. 

For Im k > the Coulomb Green's function G~^(k) turns out to be diagonal in a 
very specific Sturmian basis, namely when the parameter b has been set equal to —ik. 
Therefore, in coordinate representation, following |l5| , we have 

v— 

but it should be clearly stated that for real k the above series diverges so that ( p4) ) is 
invalid for physical momenta. Formula ( p4[ ) is quite remarkable in that the Green's 
function is given as a sum of the poles at the Coulomb bound states yet the same 
expression holds for Coulomb repulsion as well as in the Coulomb-free case. Unlike 
the function ge(k), the Green's function is analytic function of the Coulomb strength 
parameter Z. Most importantly, however, the residues in ( |24| ) factorize what results 
in a separable representation very useful in practical applications. In particular, this 
implies that the kernel of the Lippmann-Schwinger equation either (Q) , or ( |l3| ) becomes 
degenerate. With such a simplification the solution of equation (|l^) is straightforward. 



Setting k = ip with real p > and inserting ( |24| ) together with (12) in (|13|), the 
solution of ( |l3|) is 



MK r) = [1 - «(*) *<(*)] fc(*> r) - £ E S fPf X ^ p \ (25) 

p^ =Q t -f i -r v -t- rj 

where 



X iv (p) = S e „(p,r)V(r)w e (k,r)dr. (26) 



(i 



In order to determine the unknown coefficients X^ u (p), as well as the K-matrix, the 
expression (pq) for wi{k, r) is inserted in ( p6| ) and in (jlj). This gives an infinite system 
of linear algebraic equations where the unknowns are the coefficients X( v together with 
Kt(k) 



oo 



v'=0 



u (tv\V\lv') 
Ovv' H 



Xtwip) +ge(k) A tv {p) K e (k) = A lv (p), (27) 



pl + l + v' + fj 

where v = 0, 1, 2, and this infinite set is supplemented by the equation 



/X y, A lv , (p) X tv , (p) 

p ^ e + i + v' + f) [9i 



- U (k) {fa\V\M - (1/2//)] K e (k) = <&|V|&> . (28) 
In (|27]) and (|28|) we have introduced the following abbreviations 

poo 

Ae,(p)= S&,(p,r)V(r)&(*,r)dr, ( 29 ) 
Jo 

{lv\V\(.v') denotes the expectation value of the potential in the Sturmian basis 

/>oo 

(lv\V\lv')= S iv (p,r)V(r)S^{p,r)dr, (30) 
Jo 

and (</>£|V|<^) is the Born term 

POO 

(<f>t\V\4>t) = Mk,r)V(r)<f> t (k,r)dr. (31) 
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The integrals ( p9| ) - (|3l|), carrying important information about the potential, should 
be regarded as input and to proceed further it is sufficient to truncate the infinite 
Sturmian basis approximating it by a finite dimensional one. When the resulting finite 
set of linear algebraic equations has been solved by any of the standard methods, one 
of the solutions contains the K-matrix. 

The scattering matrix Tg(k) can be calculated in a similar way. By inserting ( p4[ ) 
in (§), we get 

u e (k,r) = Mk,r) - - L e+l + v + fj ' (32) 
where the new expansion coefficients defined as 

p oo 

Yt„(p)= S^(p,r)V(r)u e (k,r)dr. (33) 
Jo 

are obtained by solving the following system of linear algebraic equations 

{lv\V\lv>, Ytv ,{p) = A tv {p). (34) 



i/'=0 



p£+l + v'+f}_ 

Using ( |32"| ) in (j|), we arrive at the ultimate expression for the scattering matrix 



n(k) = -2/i 



k it/u \ ^ ^ A e ^(p)Y ev ,(p) 
<>e\V\(pe) 



p I + 1 + v 1 + fj 



(35) 



Similarly as before, for practical purposes the infinite Sturmian basis must be 
approximated by a finite dimensional one. The values of k for which the determinant 
of the system (|34|) vanishes correspond to the bound states. For attractive Coulomb 
potential (Z < 0) they are infinite in number, all shifted from the pure Coulomb 
position owing to the nuclear interaction. For Z > the bound states may, or may 
not exist, depending upon the nature of the nuclear potential. Since the case of 
bound states has been discussed elsewhere we assume in the following that for 
the considered k value the system of equations (|34|) is non-singular. Once r^(fc) has 
been determined, the effective range function may be obtained by making use of Heitlcr 
formula but this approach is equivalent to solving the system of equations ( f27[) and 
(^q). The determinant of that system may also vanish at some values of k and they 
correspond to poles of the K-matrix. Clearly, the poles of the scattering matrix are 
distinct from the poles of the K-matrix. 

The proposed scheme has, of course, its limitations and breaks down not only for 
large p but also for very small p. The large p limit is rather obvious and stems from the 
fact that for large r the regular Coulomb wave function <fie(k,r) increases like ~ e pr 
and therefore the Born term integral (|3l]) diverges for p > A/2. Consequently, the 
largest admitted p value is determined by the range of the nuclear potential. There is 
also a lower bound for p, as in the limit of small p the convergence rate of the Sturmian 
expansion gradually deteriorates eventually rendering the whole scheme impractical. 
The lower limit for p, according to our experience, depends upon the characteristics 
of the nuclear potential. In this situation, the effective range function must be 
extrapolated to zero. If the extrapolation is not satisfactory the K-matrix at k = can 
be obtained by solving the Lippmann-Schwinger equation with the Green's function 
given by (|T5|). Although the kernel is not degenerate in that case but the resulting 
integral equation can be easily solved by standard methods. Another alternative is a 
numerical integration of the appropriate Schrodinger equation. A stepwise integration 
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of the latter from r = up to some large r — R, well beyond the range of the nuclear 
potential, yields the regular solution u e (0, R) and the derivative u' e (0, R). With these 
two functions in hand, the effective range function is 

1 W[u e (0,R),8e(0,R)} 
Kt(0) W[u e (0,R),M0,R)Y 
where W[f, g] = f g' ' — g f denotes the Wronskian determinant and the Coulomb wave 
functions are given by ( |l6| ) and (|l7|), respectively. 

4. Applications 

In order to examine the performance of the calculational scheme presented in the 
preceding section, the Sturmian expansion method will be subjected to a rather 
demanding test. To achieve this goal, for the nuclear potential V(r) we have selected 
Reid nucleon-nucleon 1 Sq potential 0]. This potential is given as a superposition of 
three Yukawa potentials, has a one-pion-exchange tail and a short range very strong 
repulsive term. Since the shape of this potential shows a good deal of structure 
comprising of attraction and repulsion, the rank N of the Sturmian approximation 
must be large enough to account for the rapid variation of the potential. In addition 
to that, the other unfavourable circumstance is that the long one-pion-exchange tail 
greatly reduces the range of the admissible momenta which must not exceed 65 McV. 
In the actual computations, we have been able to determine the effective range function 
in the interval 20 MeV < p < 65 MeV. The rank of the Sturmian approximation N at 
which convergence has been attained varied from about hundred to several hundred 
at the low-energy end. It is perhaps worth noting that when V(r) has a Yukawa shape 
then all the necessary integrals are given in an analytic form. To be more explicit, 
At,, is obtained from 



Stv{o, r) <pt(k, r) dr - ' 



2 b 



{X + b) 2 + k 



Y +1 f\- 


6H 


-ik\ 


\ U+ 




-ik) 



CTL-iJ 2F 1 (-v,l+l-ir 1 ;2£ + 2;l-z) (36) 
A + o + iky 

where z = [A 2 — (b + ifc) 2 ]/[A 2 — (b — ifc) 2 ]. The expectation value of the potential in 
the Sturmian basis is obtained from 



Siv (b, r) Si v > (b, r) dr 



1 (2£+l + v)\ {2l+l + v')\ 



(2€ + 1)! V v\v'\ 
x x 2l+2 (x + iyM-2-u- v ' ^(-u, -u'; 21 + 2; x 2 ), (37) 
where x = 2b/ X. Finally, the Born term, takes the form 

<pt{k,r) 0^(fc,r)dr- 



(A 2 + Ak 2 Y +1 \\ + 2\k / 
x 2 F 1 (e + l-irj,£ + l+iT];2£ + 2;y), (38) 

where y = 4fc 2 /(A 2 + 4k 2 ). In the above expressions we assumed that A 2 > 4fc 2 and 
2-Fi(a, b; c; z) denotes the hypergeometric function. 

The effective range functions resulting from our computations are presented in 
Figure 1 as functions of the cm. energy E — k 2 /2/i. The calculations have been 
carried out for three values of the Coulomb strength parameter Z: Z = ±1 and 
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Z = 0. Strictly speaking, attractive Coulomb potential would be appropriate for 
the proton-antiproton pair. However, in that case the Reid potential, devised to fit 
the nucleon-nucleon scattering data, is not applicable and the inferred effective range 
parameters have no relevance to proton-antiproton scattering. Nevertheless, this case 
is of theoretical interest and has been included here for testing purposes. In Figure 1 
the continuous curves above threshold result from a numerical integration of the 
appropriate Schrodinger equation. The entries below threshold, depicted by asterisk 
symbols on the plot, have been obtained by using the Sturmian expansion method. 
As expected, close to E = the values of the effective range function calculated below 
and above threshold appear to be lying on the same straight line. In order to qualify 
that statement, the values calculated below threshold have been subsequently used 
to extrapolate each of the three functions to zero. This gave the intersection points 
at zero together with the slopes of the curves. The scattering lengths and effective 
ranges, extracted that way with the aid of Jig) , are presented in Table I (referred to as 
"below"). For comparison, we give also their counterparts calculated from the above 
threshold solutions of the Schodinger equation (referred to as "above" in Table I). The 
agreement between the below threshold and the above threshold computations is in 
all three cases very good. The error on ro is slightly bigger than that on ao, as the 
computation of the former involves the differences of the function values what results 
in some loss of accuracy. The effective range seems to be little affected by the Coulomb 
interaction what is also reflected in Figure 1 in that the three curves have almost the 
same slope near zero. By contrast, the Coulomb corrections to the scattering length 
are quite large especially for Coulomb attraction when ao has been pushed far away 
to positive values. The K-matrix, unlike the effective range function, shows rapid 
variation as function of energy and in all three cases exhibits a pole not far from 
the origin. At the pole the effective range function goes to zero, so that the actual 
positions of the poles of the K-matrix are apparent from Figure 1. It is interesting 
to note that for attractive Coulomb potential the pole of the K-matrix lies in the 
physical region, very close to threshold. This example is rather instructive, because 
it shows that any phenomenological attempt to approximate here the K-matrix by a 
polynomial in k 2 , with coefficients to be fitted to the low-energy data, would have a 
disastrous effect. The problem is immediately alleviated when instead of the K-matrix 
the effective range function is employed. 

The two-potential problem considered in this work can be easily generalized to 
comprise complex nuclear potentials. Such potentials appear in various problems of 
nuclear physics and, in particular, they have been widely used for a phenomenological 
description of the hadronic atoms phenomena. The introduction of a complex nuclear 
potential into the description of a hadronic atom has two observable effects: a shift 
e in the energy of the otherwise hydrogen-like level position, and a broadening T of 
the atomic level arising from the fact that the system can now decay through the 
absorptive strong interaction. When the orbital momentum is large enough the level 
shift in a hadronic atom is usually quite small in comparison with the Coulomb energy. 
This observation has prompted several authors (cf. (l^, |l9[ |l(| [2l], |2^]) to introduce 
small shift approximation (SSA) based on the idea that if atomic bound state pole 
has been shifted from the Coulomb position fc„ = — i fia Z/n, where n is the principal 
quantum number, by 6k then \Sk/k n \ can be regarded as a small expansion parameter. 
One of the implications of SSA is the shift formula expressing the energy shift in a 
hadronic atom in terms of hadron-nucleus low-energy scattering parameters. This 
formula relating only experimentally measurable quantities without reference to the 
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underlying nuclear potential has very important phenomenological applications. 

In order to derive the shift formula we are going to obtain an approximate solution 
of the bound state equation (|l9|). As mentioned before, the function ge(k) that enters 
( fl9| ) is singular at k = k n and for k such that \k — k n \ is small in comparison with 
the separation of the unperturbed levels \k n — k n +i\, we can approximate ge(k) by a 
single pole term 

«(*) « ■ (39) 

ft ft n 

The residue i?^ n can be worked out from the definition of gi(k) and reads 

Ren = { 2 \ k nf ^(fc„,r) 2 dr} * (40) 

where the normalization integral, is 

, s 2 , n(n- £-1)1 [{2£+iy.\e\} 2 

**k»,r)**r= 4|fcra|2 , +3( , + n) , • (41) 

In order to obtain an approximate solution of equation (|l9| ) it is assumed that the 
left hand side is a slowly varying function of k so that 1/Ki{k) may be replaced by 
1/Ki(k n ). Using the single pole approximation (|39|), the solution of (|l9| ) is 

k^k n [l + K e (k n )R en ], (42) 

and this yields the level shift formula in its most general form 

8E = -e - i|r = E - E n k -K e (k n ) {2^ /" 0,(fc„, r) 2 dr}" 1 (43) 

where E n = k n /2/j, and it is understood that for complex V(r) the K-matrix is 
necessarily also complex. Approximating in ( (43| ) the K-matrix by the scattering 
length one obtains the celebrated Deser-Trueman formula jl|, ^9|. However, the latter 
approximation is hardly necessary as the K-matrix below threshold can be calculated 
with the same amount of labour by using the Sturmian expansion method described 
in the preceding section. The other alternative, that does not require the SSA ansatz, 
is to locate the poles of the scattering matrix. Using Sturmian expansion, the problem 
reduces to finding the values of k at which the determinant of the system (|3~4| ) vanishes 
(cf. J!!). 

In conclusion, we have presented a simple and efficient method for calculating 
the effective range function below threshold and the same calculational scheme can 
be used for solving the bound state problem. The main advantage of the proposed 
method lies in the exact treatment of the Coulomb interaction what is accomplished by 
using the Lippmann-Schwinger equation whose kernel involves the Coulomb Green's 
function. Since this kernel below threshold assumes degenerate form, the solution 
of the integral equation is straightforward. The essential point is that the Coulomb 
Green's function and the regular Coulomb wave function, constituting input to the 
Lippmann-Schwinger equation, are real analytic functions of k 2 , free from Coulomb 
singularities. Both are given in analytic form suitable for numerical computations and 
the irregular Coulomb wave function is never needed. 
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Table 1. Coulomb-modified low-energy scattering parameters for Rcid 1 So NN 
potential calculated below and above threshold. All entries are in fm. 
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Figure 1. Coulomb-modified effective range function vs. cm. energy for 
Reid 1 So NN potential. The different curves correspond to Coulomb attraction 
(Z = — 1), Coulomb repulsion (Z = 1) and Coulomb-free case (Z = 0). 



